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XXVI. — On the Rotatoby Motion of the Heavenly Bodies.- By 
the Bev. W. G. Penny, M. A., Professor of Mathematics in the 
Catholic University of Ireland, and late Mathematical Scholar in the 
University of Oxford. 

[Read February 24, 1868.] 

1. The object of the present Paper is, in the first place, to ascertain 
whether the various disturbing forces which act upon the heavenly bo- 
dies produce any permanent effect upon their rotation ; and, secondly, 
supposing such an effect to exist in general, to ascertain under what 
circumstances it will cease to do so — that is to say, what are the con- 
ditions under which they would rotate permanently without any but 
periodic changes. 

The inquiry may be of some interest, for two reasons : — First, it has 
been ascertained that the observed acceleration in the motion of the 
moon has been only partially accounted for by the diminution of the 
excentricity of the earth's orbit, to which cause part of it, though not 
much more than half, is undoubtedly due ; that is, if we calculate 
what onght to have been the angular distance of the sun and moon at 
the time of an ancient eclipse — say 2500years ago — it is found that, 
after making allowance for the acceleration produced by the cause men- 
tioned, that their angular distance so calculated does not agree with 
what it was actually observed to have been. Now, such an error might 
be produced either by an error in the supposed velocity of the moon, or 
of that of the sun, i. e. of the earth ; but both these have been care- 
fully examined, and found to be inadequate to explain the phenomena. 
There is, however, a third cause which would give rise to the same dis- 
crepancy between theory and observation — namely, an error in the mea- 
sure of time — just as an error in a ship's longitude at sea might be 
caused by an error in the rate of the chronometer. Just so an error in 
the angular distance of the sun and moon a certain number of years ago 
might be caused by an error of any one of the three elements which 
enter into it — the length of the year, the length of the month, and the 
length of the day; and so it is evident that, if the length of the day has 
undergone any sensible alteration since the date of the eclipse above 
spoken of, the real length of time that has elapsed since then must be 
different to what would be supposed if the length of the day had re- 
mained invariable ; and thus the relative positions of the sun and moon 
would also be different. Accordingly, it was suggested by M. Delaunay 
that possibly the length of the day — or, in other words, the velocity of 
the earth's rotation — had varied. 

In the next place, the moon, and it is said satellites generally, turn 
always the same face towards their primaries ; that is, they rotate about 
their axes in the same time that they perform a revolution about 
their primary. And it has been asked whether there is any cause for 
this ; in other words, whether the body about which they revolve exerts 
any influence upon them which would affect the rotatory motion, so as 

a. i. a. peoc. — vol. x. 2d 
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to make it coincide in period with that in the orbit, supposing that 
there had ever been a time when they did not do so. 

The general results arrived at may be thus stated : — 

(1). That in all bodies which are perfectly symmetrical with respect 
to the three planes containing their principal axes of rotation, such as 
an ellipsoid (not one of revolution), there is no permanent change pro- 
duced; nor in any in which the moments of inertia about the two 
principal axes which are perpendicular to that about which it is revolv- 
ing are equal; I think there is a permanent change whenever the three 
moments are unequal, though the body may be perfectly symmetrical 
with respect to each axis ; 

(2). That in other cases there is in general a permanent change 
produced ; 

(3). That in the case of the earth disturbed by the lunar moon, the 
change produced will be so very small as to account for a very minute 
fraction of the whole amount required to explain the phenomena above 
alluded to ; 

(4). That one condition under which there will be no permanent 
variation, is when the time of rotation nearly or exactly coincides with 
that in the orbit ; but this is only one out of several other such rela- 
tions as might exist; just as there are always several positions in 
which a body might remain in statical equilibrium ; and that in some 
cases, though not in all, the forces are such as to produce the relation 
above spoken of; and, lastly, 

(5). That the effects are enormously more rapid in the case of a sa- 
tellite described by its primary than vice vend. 

I have supposed, in treating the question, only one disturbing force 
to be acting upon the body, and also that its orbit is a fixed plane ; 
neither of which, especially for the earth or moon, is strictly the case, 
but will be sufficiently near the truth for the present purpose ; also, I 
have supposed the body to be entirely solid, instead of being partially 
covered with a thin layer of fluid. Mr. Airy, however, is said to have 
examined the effect of the tidal wave, which it was supposed might by 
its position, &c, produce some retardation upon the earth's motion, and 
has found it to be insensible. 

Eesult (4.) has been spoken of as only an approximate one — indeed, 
to pretend to extract anything more out of differential equations which 
can only be solved by successive approximations, as is the case in the pre- 
sent instance, would seem almost to amount to a contradiction of terms. 

2. The differential equations of motion are 

dw l C-B 1 *yi - yz i ■» 



dw, A-C 1 „ x*\- »i 

+ - IUjO). = — 2.11% = "= 

dt B B [x-xf + y-yf + Z- 

dw, B-A 1 yxx -*y. 



A. 
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The solution of these equations will give the values of to ,, u> 2) and oj 3 , 
and these being known, the value of w = */us? + ">.' + °>*> or the value 
of the velocity about the actual or instantaneous axes of rotation is 
known also ; and if any one of the quantities a> u <w 2 , or io 3) contain any 
term which increases with the time, and is not periodic, there will be 
a permanent change. The body is supposed to be rotating very nearly 
about the axis of z, about which <o 3 is the velocity, and where (■>„ a> 2 , <u 3 , 
are the angular velocities about the principal axes, A, B, C, the moments 
of inertia about the same; also x lt y t , z u are the co-ordinates of any par- 
ticle m referred to the principal axes, the origin being at the centre of 
gravity of the disturbed body, and x, y, z, are the co-ordinates of the 
disturbed body, the plane of x, y being that of the orbit, the intersec- 
tion of the planes of x u y„ and x, y, the axis of x, and the origin as 
before. 

Let i be the angle between the planes of x, y, and x u y, ; the lon- 
gitude of the disturbed body measured from the axis of x on the plane 
of the orbit, the right ascension of the axis of x l measured on the 
plane of x lt y, that is, the angular distance of the axis of x x from that 
of x ; r the distance of the disturbing body from the centre of the dis- 
turbed, and Ti the distance of any particle m from the origin. Then we 
shall have by spherical trigonometry, 

X 1 + cos » . , „ 1 - cos « , , „. 

- = COS (0 - 0) + — cos (<p + 0) 

r 2 ■* 

y 1+cos* . ., _ 1-cosi . 

*= — sin(0-<?) r — sin (0 + 0) 

r £ * 

% = - sin i sin 0. 
Equations which are usually given in the form 

- = sin cos + cos i cos 9 sin 0, &c. 
r 

But the form given above will be much the most convenient for the 
present purpose. The equations (A) can only be solved by successive 
approximation. The first approximation will be when the right hand 
member is 0, that is, when there is no disturbing force, or when the dis- 
turbed body is spherical. The next will be when the bodies are supposed 
to be spheroids of revolution. This very nearly represents the case of 
the heavenly bodies; but inasmuch as they have a variety of irregu- 
larities both of form and density, will not accurately do so ; and it be- 
comes therefore necessary to examine what will be the general effect of 
the said inequalities of surface ; and specially to see whether there will 
be any permanent alteration in the velocity of rotation arising from 
them. Supposing such to exist, it is manifest that in consequence of 
the bodies being so nearly spherical, it will take place very slowly ; 
but the ultimate amount of alteration will be none the less than if the 
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inequalities which produced it were more considerable ; only it will 
take more time to arrive at a permanent state. 

To calculate, then, the effect of the disturbing force, we must de- 
velop the right-hand members of equations (A). To begin with the 
first. The quantity to be developed may be put in the form 

2m (ay, - ya,) { (a* + y* + a') - 2 {xx x + yy x + aa,) + ar, 1 + y* + *,*}-* 
or _ Sm (ay, - ya,) ^1 - 2 + -J 

= - 2m (qr, - ya,) ^1 + 3 - ^ + 

1 g 

= _-2>» (ay, - yz,) + -r 2m (sy, - ys,) (ax, + yy x + 2a,) 

1 3 

+ ay 2m (y, ! - z, 2 ) + (s* - y 2 ) 2m ( y,s, - *y 2(m 2,3:,) 

rejecting for the present the further terms in the development. 

By the property of the centre of gravity and the principal axes, the 
terms 2(m y,), &c, and also 2(m a^y,), &c, vanish ; and it is re- 
duced to 

3 3 

— zy 2m (y, s - a, 5 ), or — ay 2m (y, ! + x* - (a,» + a;,')) 

or ^ay(5-(7), 

hence the first of equations (A) becomes 

dw. B-B 3fi C-B 

: + OMB, = ■ %V ; 

dt A 3 2 r* J A 
similarly dw a A - C Sp A - C 

— — + =; — w.w, =— %X — 

dt B ' l r> B 

dw 3 B-A 3/t B-A 

and if in these we substitute the values of xy and a, given above, they 
will become 

dw, C-B % ix C-B . / 1 + cos « 

a } + asiu, = - sin « — cos (<p - 26) 

dt + A ' 2 2 r 3 A \ 2 ^ ' 

1 - COS « . 

cos / cos — cos (<p + 20) 
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dw a A - C 3 /i A- C . fcoa i + 1 . 



<?w 3 2? - ^ Z ft B- A 

dt * C ' 2 2H 



. , cost - 1 . , , „„) 
cos » sin + — sin(0 + 20) j 

- A il + cos t 2 . ._ , _„ 
^- | 2 B1I1 (20-20) 

1 — COB ( ' 1 

+ £ sin *< sin 20 + sin (20 + 20J 



In the first of these three there are three terms on the right hand side, 
each of which being integrated will give a term in w, of the general 
form H sin (x) ; and corresponding to this there will also be a term in o>, 
of the form JTcos (x)- Now, if these two are substituted for ui, and 

B - A 

«Oj in the function — ^— u>\u> % which occurs on the left side of the 


third of the above equations, it is evident that they will only produce 
a periodic term ; but this is because one of them is a sine, and the 
other a cosine. If, however, they had both of them been sines, or 
both cosines, the case would have been very different, and the multi- 
plication of them together would have produced a constant term. 
What we have to do, therefore, is to see whether the further develop- 
ment of the disturbing function will produce any such terms. And it 
is very readily that it does produce a considerable number of them, cor- 
responding to different combinations of and of the Bame kind as 
those in the equations last formed. Those which I shall select for exa- 
mination at present are those which have <}> - 6 for their argument ; so 
that we must develop the disturbing function so as to include all terms 
of the form sin - or cos 0-0 wherever they occur in the first tw° 
equations ; that is, those for <u, and w % . 

Let us resume, therefore, the two last terms in equation (0) which 
had been rejected, and which will contain all the terms of lowest di- 
mensions of the form required, we shall then have for this part of the 
function on the right side of the equation for <t>, 

1 p. ( 3 „ , .r, s 15 .(xx, + yy, + as,)* 
A ? | 2 ( * y ' " ySl) ?" T (8y ' " V ^ ? 

substituting the val ues of z, y , &e., in this, and retaining only terms of 
the form sin0- 0, cos - 6, the first part of this expression will give 



for shortness, let 



13ii v , , . 1 + cos • . . 



1 + COS 1 1 - COS J 

- — = «and — =/J 
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also let 2 (»»>,'«,), &c, be denoted by rfa, and similar expressions for 

the others. 

Also the latter part will give 

- — — - 2« (ay! - yzi) (izV + y'y> + a»»,« + 2xy x x y x + 2xz x x z x + 2y*y,s,) 

./l it T 

Multiplying these two factors together, and retaining all such terms as 
are either of one dimension x and y, or of three dimensions, which are 
the only ones which will produce terms of the form required, we 
have 

- -j — - ( afy, xfr + (y 3 - 2«V) y?*i + #y *i 3 + H*f - x * % ) x $i z \ 

J&. Jt T \ 

Now, 
x = a cos - + p cos + 0, and y = - a sin - - j8 sin + 6 

.-. x 2 = ± (a'+J/?) + £a* COS 20- 10 + a/9 cos 2?) + a/9 COS 20 

+ \§P cos 20 + 20 

Multiplying this by the value of y, and retaining terms whose argu- 
ment is - 0, 

^V " (~ i( flS + a/ 3 *) + 4°* + aW - i / 3 *) sin 0^ = - ±(a 3 + 2/9*a) sin 0^0 

which is the coefficient of x x *z x above. 
Also in like manner 

y» = - 4 (a 3 + 2/3 l a) sin (0 - 

and 

sin' i ** —— — *— ^— 

* 3 y = — - (1 - cos 20) (a. sin - + sin + 0) 

= - i sin'< (2a - /9) sin - ft 

Thus the coefficient of y?z x becomes 

- (| (a s + 2/8 8 «) - sin'. ) 4a - 2/9)) 

and that of 

ij 3 is - | 6in*« 2a - /9 sin - 0, 
also 

y* = £(a 2 + /9 s ) - ia s cos 20-20 + £a/9 cos 20 - £a/S cos 20 

- J/8 2 cos 20 +"20 

.-. 2y* = A(a 3 + 2a/3 ! ) COS 0-0 
and 

rz 2 = ^ sin s i 2a - /3 cos - 
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so that the coefficient of 

a^z, is %(a 3 + 2a/3» - 2a - /3) cos <j) - 
making these substitutions, we have 

+ J T V r + 2/3,fl ) *■**> + W" 3 + 2/3 *°) ~ 4a ~ 2/3 sin,t ) ^ 

+ sin** (2a - 0) *?} sin (0 - 0) 
1 15 fi 



, — , (a s + 2(1/3* -2a- p sin*«) xtf a cos - 
^4 4 r 

and the whole coefficient of sin - will be the quantity just found, 
minus the quantity first found above, that is, it will be 

J- 2 \y I (° 3 + W ^ + ( 3 (° 3 + W sin2 « ( 4a - 2Z 3 )) yi^ 

3 ) 

+ sin»« (2a - 0) z, 3 j - - a r^, 

call this, for shortness, ~j My, and the coefficient of cos <p- 0, just 

given, - — N. Then the equation for a>, will become 

rfo., C-B 1 „, . -j—, 1 , T -r— » 

— + — - — jojioj = -j-M, sin - (? - T iv cos 0-0, 

and in like manner 
dw t A-C 



,-- + — 5— a>jo», = - — N sin <b -0 + - i rM x cos<f>-0 
at a a is 

where M, is what M y becomes when x x and y, are interchanged. In 
order that these may be integrated, it will be necessary to express and 
in terms of t ; but before doing so it may be well to make a remark 
upon the quantities 2(i» xfa), &c, which occur continually. We may 
always choose the axes so as to make 

2 (mtd) = 0, 2 (mxtfj = 0, &c. ; 

but not so as to satisfy any further conditions, such as 

2 (mx.%), 2 (mxfa = 0. 

But there is one case in which a certain class of these quantities will 
always vanish. Whenever the body is perfectly symmetrical in term 
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and density with respect to the planes which contain its principal axes, 
it is manifest that all quantities of the form 2 (m *,%«), &c, when 
either p or q are odd numbers, will vanish. Thus, suppose q to be an 
odd number, then any particle m having s, for one of its co-ordinates 
will always be accompanied by another having - 2, for one of its co- 
ordinates; hence the sum of these will vanish. Now, these are the 
sort of quantities, both in the example which I have chosen, namely, 
terms depending upon sin - 9 and cos - 0, and in all other cases 
whatever where constant terms can be produced. Where this is the 
case, therefore, the quantities in the equations for u> lt &c, on the right 
side will vanish ; but this will not be the case for Buch bodies as the 
earth, moon, &c, whose form, though nearly spherical, &c, differs 
from it on account of the irregularities of surface, &c. 



Integration of the Differential Equations. 

For this purpose the following equations for determining tp, &c, 
which are given in all dynamical treatises of motion about a fixed point, 
will be useful, viz. : — 

d<j> cos« . 

— - = « 3 — : — (», sin <p + «j cos <p) 

at smi 

tff _ 1 



dV 1 

—L- = — . — f„ s in <s + « 2 cos p) 

at sin « 



di 

-j — «! cos <p — aij sin <p, 

in which v[» is the longitude (measured backwards) of the moveable axis 
of x, referred to a fixed line. 

The first of them will give us an approximate value for : to use 
it we must first find «e 3 . Now, neglecting the disturbing force, as also 
products of <Ui«u 2 , which are supposed very small, and are, moreover, 
multiplied by B - A, the third of equations (A) becomes 

dot- 

— 3 = 0, or » 3 = n 

Neglecting, therefore, small quantities in the first of the equations of 
this number, it becomes 

— = n, or <p = nt, 

if the time is supposed to commence when the R.A of the example x is 
; also for 6 it will suffice to put », + «. Thus the equations for *>, 

and *>j become, on putting n for « 3 , nt for <p, &c, 



197 
-p + — ^— m, = -rM s sin (»"», t+e)- -r N cos (» - «! t + «) 

*!» + £^ ws , 2 =_I iVsin(^¥ l *+ e)+ I Jf, COS(?T^; < + e) 

ill li Ji -D 

If we differentiate the first of these, and for the value of — 2 , which -will 

occur in the result, substitute its value as derived from the second, we 
shall obtain the following equation, from which *, has been eliminated: 

n-rf C - Bn\ . 



+ N\ —j— + ~. B - sin (n-n,t + e) 
The terms introduced into m, by the integration of this, will be 
M , G - B nr 



c-bc-a ; ; — ; C08( - n ~ nH+e) 



,,,.-», C-B 
N[ -J- + ^B> 



- C-BC-A , f -"»(-".« + .) 

— _ — »« - (n - »')* 

call this, for shortness, 

At cos (»-«!<+ e) + Bi sin (n - n, tf + e) 
In like manner the integration of the equation u> 2 will give 

«>, = .4 2 cos (»-n, t + e) + Bt sin (re - wi t + e) 
where 

A * = ~ C-BC-A . . B *"~ C-BC-A 



The complete integral will also contain arbitrary quantities of the 
form 



[C-B C-A \ 



1. I. A. PEOC. VOL. X. 



6 COS 

V \ A Ji I 

2i 
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For the earth, however, the constant e is quite insensible ; but for the 
moon, Laplace says that it is variable. However, it will not affect the 
recent inquiry in either case, and therefore may be dismissed for the 
present. 

Formation of the Constant Quantities in the Differential Equations. 

If, now, we multiply together the two values of w l and <uj, given 
above, we shall have 



»,«, = [A x cos (m - », t + e) + B x sin (« - n x t + «)) 



{^2 cos (n - »i t + e) + B 2 sin (n - n x t + e) j 
= \A X A^ + B&) + 
periodic terms ; retaining the constant part 
B-A 1 B-A, 



•,i, = - 



2 C 



A X A % + B X B. 



or, putting for A u &c, their values, and dividing numerator and deno- 
minator by 

— — , C-AC-B 

n ~ n AB — 

1 B -AN.(M y +M x ) 



where 



2 45C D. 



_ C-A C-B 2 . 

D = 7-= n 2 - n - »,' 



This part, therefore, of the differential equation for <u 3 contains' a 
constant term ; but before we can say that the entire equation does so, 
it is necessary to develop the term on the right aide of the equation. 
Now, if this is expanded, it will be easily seen to consist entirely of 
sines and cosines, of which the general form may be said to be 



/ sin (pn - qn x t + a-) 

where I is some function of i, and p and q are whole numbers. It 
would appear, therefore, at first sight, to contain no constant term ; 
but in reality it will be seen that it does. For, it is easily seen that 
every term in <u, and a> 2 , such as those found above, will introduce a 
periodic term into the value of «, as also of $ and 8, and the multipli- 
cation together of periodic terms may produce a constant To see what 
terms in the development spoken of will be necessary, we must find the 
variations of *', <{> and 0, then, 
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di 

— = «r, cos ffl - s>j sin 0. 

at 

Put for its first approximate value nt 



*i cos = •§ (^4, cos (2» - re, < + e) + 5, sin (2« -«,« + «) 

+ -4, cos (»,£ + e) - J5, sin (nt + *)) 



» 8 sin<p =i(4 2 sin(2n-»,« + «) - # s cos(2n-n 1 < + e) 

+ -4 a Bin (»,£ + «,) + i?* cos (»,< + «)} 

taking the difference, and integrating, 

l = { +l( J h±Jh sin(2^«;« + e) + ^— - COS(2M-n,< + e) ] 

• ¥ \2»-Mi 2»-», y 

+ i _! ? sin OM + «) + c° 8 (M + € )1 

8 V »i n i ) 

for the variations of we have 

rfffl cost . . . 

-J- = », 3 - (ar, sm « + »j cos a), 

at sin ( 

its variations, therefore, will arise partly from those of v„ and partly 
from those of a, and » v "We may set aside the former for the present, 
and confine ourselves to those of «, and « 2 . 
Then, priming the function 

cos « , . 

-; — (*>i sin + «j cos ffl) 
sin < 

exactly in the same way as above, and integrating, we have 

, Icosj {At + Bi At- By . ) 

a = nt + {—- cos (2» -»,< + £ — Bin (2w- 71, J + e)\ 

f 2sint(2»-», v 2»-», v 7 J 

lcos« L4i--8 2 , . . A'+By A 

+ « -r— 1 cos (n t t + e) ; — Bin («,< + e)J 

2 sin i ( re, v i y w , /j 

also since the differential equation for ^ is 

-i- = - - — (*>, sin a + #j cos p) 
al sin * 

the variations of ty will be the same as those of 0, only not multipled 
by cos i. 
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We must, therefore, carefully seek out all terms in the development 
of it, whose argument is either 2<p - or 0. The latter, however, will 
be found to disappear. Then, 

= &~ 2 ^ ~ ^ ~2~ 1 CC * a;,S + 2/Sy ' 2 + **** + 2Xy ^ + 2a * *'*' + 2y,? ^ A 
= -p- 15 ^ 'W(y ! - *> + (a;, 2 - y. 2 ) ay*j 

= + — -^- 15 ( o s cos 2p-20 + 2a/3 cos 2p + /3 s cos 2p + 25 



x yfa sin « sin 



1 /•!«/' ,.:. 



^ _ . fL — I „? sin 2<p - 20 + ap sin 2<p + /S 2 sin 2$ + 20 {x*-y*) 
u r^ 2 ^ y 

sin < sin 

= — — — (o 2 - 2a/3) a;,y.^i sin « sin 2p — 
C r* 2 



sin i = sin 



- — 4 -T- (<"•' - 2a/3) (a:,* - y, 2 ) sin i cos 2<p - 0) 
r 4 

Now, taking the value of i found above, we shall have 

X A^B 2 . \ 

^ + ^2^^ Sm2 «- w ' + -) 

neglecting, for the present, the other term, 

, A l + B s . 

= sin «, + + cos «, — sin 2re - », + 

* 2n-n l 

A + B 2 . 

cos i = cos «i - -k sin t, — sin 2n - n, -, &c. 

* 2n-n 

Then, since a stands for — , we shall have 

a? = i f 1 + cos *! 2 - 1 + cos «, sin i, — r 2 A, + B* sin 2» - n 1 ) ) 

\ 2n-n' J) 



and 



A + B. 
- 2afi = - i sin 2 « = - | ( sin 2 <; + sin d cos «, — * sin 2% - n } 
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.-. a s - 2a/3 = i (1 + cos «« - 2 Bin 2 *,) - J sin «, (I + 3 cos «,) 

2»-»' 
Multiplying this by the value of Bin *', we have 



sin 2n - »', 



(a* - 2a/3) sin i = £ (1 + cos i,* - 2 sin**,) sin <, 



+ 1 1(1 + cos i,* - 2 sin ! 4,) cos «, - £ sin *»,' (1 + 3 cos «,) 



A, + _B 2 . 



, sin2n-»,. 
In - n' 

Also, considering the variations introduced into <p, we have 

, cos «, ^4, + -B s 

r * sintj 2»-n' 

also is the true longitude measured from the moveable axis whose 
longitude measured backwards is yr ; therefore if »,£ be the longitude 
measured from a fixed axis, that measured from the moveable axis will 
be n l t + V" ; or, from the value of ty given above, 

u 1 1 A 1 +B 1 

nH + - - r cos 2» - w, 

2 sin *, In - n l 

2 008^-1 A t + Bt 



.:2$-0=2nt-nH + i : T coa2n-n\ 

r * sine, 2n-w 

and 

. z , , 2 cos «, - 1 A. + B» 

sin2ffl-e = sin2»-w 1 + i ^ -! A. 

Multiplying thiB by the value of (a 3 - 2 a/3) sin «, we shall have for the 
constant term depending upon A t + 2? s , 

jJL ( l+cos^ 2 - 2 sin ».,) (3 cos «, - 1) - ^ sin \ (1+3 cos i.) j ~~r 

which may be put into the simpler form 

- 1 + cos », -- sin*<, (1 + 3 cos «,) — 

4 \ 4 J 2n-n l 

and in like manner the quantity (a 8 + 2a/3) sin « will contain the con- 
stant quantity 



K 



5 . , \ A - A 

1 + cos i, - - sin**, (1+3 cos O -r 

4 v ' } 2»-» l 
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On substituting these values, the right-hand side of the equation for 
determining » 3 contains the constant terms 

_^_^l + co Sl --sin^( 1+ 3cos«)j^,-2^r^r 
This latter factor may be more conveniently put into the form 

Let us first find the value of the first pair of these terms. Now, 

and 

_ 1 „/»-«' C-J5 



where 
Jkf„ = 



il |i£ ( („» + 2/3 2 a) «,*«, + 3 (a 3 + 2/8*a) - sin s i 4a - 2/3 j y,% 



+ 3in 2 « 2a - /3 z, 3 - -= ar.'zi 



2V = — — (a 3 + 2a/3* - 2a - /S sin 8 <) *■$& 
4 r* 

then , 



becomes, considering first the term multiplied by n - «', 

jj!L. |1£ f ( a » + 2/3*0)"^, + (3(a 5 + 2/3»a) - sin n 4a -2/3) ^ 

+ sin 2 . 2a-/3 Zl 3 j -jirSiiifA+^W' ,i+ ^ a 

. ) n — n l 

- la - p sin *<) (*,»*• - y,*«,) -^p- 

or 

ii. ( 1£ | (2(a 5 + 2/3 J a) - 2a~^"/j sin h) (i^i, +#¥,) + sin *c 2a - 0*. 

3 —j-X n-w 
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This coefficient, it is easily seen, is equal to £ {M, + M y ), so that we 
may replace the above expression by 



% M x t M s a;,y,s, 



n- n> 



2 jo-z T -Uly -,,,_, 

in like manner it will be seen that the part multiplied by n becomes 

1 M s + 3t. B-C 

2~D- X ^-AB n 

thus, the two terms in question become 

\M X + M yl #,y,z, (n-n, B-C 



2 D 2n-», V A AB 

and in like manner the latter pair become 

l(M x + M y \ x,y,z, I n-n} A-C 
"2~V B ) 2~n^n~, [~B~ + ~AB " 

therefore the sum of the four is 

1 M, + M y x,!/^ f- — -/ 1 _1\ _ B-B 



D 2» - »n \AB AB 



n - n\ — - r=r 4 — — — n 



IB-^A 
2 ^5 



M x + M y x x y^ 



tnus the constant terms of the right-hand side of the equation for a 3 
become 



— 1 + cos i - — sin 5 i 1 + 3 cos i'-" * ' — 



■■) *i«/i«i 



r> 8 V 4 ; "" ' 2 ^£C V D 

also on putting for a and /3 their values, it is easily seen that the value 
of iVis 

15 fi 



— — I 1 + cos i — sin h 1 + 3 cos * ^lS^i 
8 r* V 4 / 

so that the expression just found for the constant term in the deve- 
lopment of N reduces itself to 



\ B-A I M„ + My \ 
2ABC\ D ' 



And this is identical with that part on the other side of the equation for 

da 

~, which arises from the multiplication t 

lit 

these terms identically destroy each other. 



du 

— ', which arises from the multiplication together of *, and * 2 , so that 

lit 
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No account, however, has as yet been taken of the first terms which 
occur in the development of L, namely, 



3 /* B-A 
2,1* C 



ja 2 sin 2<p - 20 + 2a/3 sin 2<p + /3 s sin 2p + 2d\ 



Now, it is evident that these will contain constant terms, which, if 
they do not identically destroy each other, or are destroyed by terms 
which may arise in other parts of the differential equations, will give 
terms indicating a gradual and permanent change of motion. Let us 
see how constant terms might arise in the above expression. 

In the first place it has been already seen that the expression for t* 
contains, amongst others, terms of the form 



ZTsin 2» - »', JTcos re 1 , &c, 

we shall thus have 

cos * = cos (t, + _H"sin 2» - n 1 + iTcos »,) = cos («, f P) 
suppose 

P 3 



cos *,{1 -i-P+7} +sin«,[P- — "Y 



= cos i t { 1 - -H"sin In - n' JTcos re 1 j + P sin «„ &c. 



= cos «, { 1 - i HKsm In - 2n l + £ HKsm 2») + P sin <!• 

In like manner if sin (2<p - 20) be developed, it will contain the 
term sin 2» - 2n\ terms multiplied by P, and constant terms. The 
first of which, when multiplied by HK sin 2» - 2» l , will produce a 
constant term, and the second when multiplied by P sin »'„ and the 
third when multiplied by cos i,. 

Again, the first approximate values already found for »,p and 0, 
will, by substitution in the differential equations, produce terms having 
the arguments 

2» - 2re„ 2m, &c, in i <p and ip-. 

And these will arise in two ways : first, it is evident that such a term 
as S cos re 1 occurring in the expression for i, would, when introduced 
into a or cos i, where it occurs in the terms in the differential equations 
already used for determining terms in «„ &c, having the argument 
tp - or n - re 1 , would introduce terms into *„ &c, having the argu- 
ment n - 2»„ and these, when multiplied by cos p, where it occurs in 
the equation for i, viz. u l cos <p, would produce terms having 2re - 2re' 
for their argument in the expression for i. These terms, however, will 
be multiplied by higher powers of sin * than those which arise in the 
manner about to be examined, and therefore for a first approximation, 
at least, may be neglected, especially in cases where »' is small ; and 
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even where it is not, it will be best to reserve them, and, if necessary, 
to take them into account after lower powers have been examined. 

The other, and more important way in which such terms arise is as 
follows : — Take the value already found for <p„ viz. — 

. 1 cosi, iAi-B-i ; A x -B t 

nt + - - . I- , cos In - a 1 + ; — cos a 1 

2 1 sin t, (2» - n* a 1 

A,-B x . - ; A 2 + B l -. ,) 

- sin 2n - n l — sin «'J 

2ii - »' a' ) 

this will give 

1 cos^ ifAi-Bi A x + B t 



cos 



, 1 COS *! ( (l 

= COS»- - — {I - 

4 1 SHli, (\ 



2»-a 1 



sin a - a 1 - r + : — 008 n - a 1 

\2a-a 1 n l 



Multiply this by »! = A x cos n - n l + B x sin n - a 1 , and retain the terms 
involving cos 2m - 2a 1 , and we have 

1 cos«, l(A,-B 2 A.+BA 

Oil cos = - . ■ { ; -x f *! 

8 lsini! (V n 2n-n l J 

{At -By. A. + B x \.\ — - 

+ ~ p + . )A l \ cos 2» - 2n' 

\2n-n l a 1 j ) 

in like manner we shall find 

. lcosi, l(Ai+B t A 1 -B 1 \ 

Ul-i BID = — { — r + ; — -"*» 

8 smt, (\2a-a 1 «' J 

A-.-Bt A 2 +B 
2a - a 1 a' 



.B 2 > cos 2» - 2a 1 



di . , lcos« (4, - i? 2 A - 2?i ^, + ^,4,+JJ, 

.-. — = a>,COS0-tt>jSin0=- ( : + 



it ' ^ ! v 8 sin i 2» - a 1 a' 



a 1 2» - a 1 



cos 2ra - 2re' 



.-.« — *,+-+ --r— - f^T+lT, ^i + j9 2 -^i--B 2 ^ 2 -^i[ 

8 sin < ( ) 

— — ; sin (2« - 2k 1 ) 



(2» - a 1 ) a 

B. I. A. PltOC. — VOL. X. 2 F 



206 

To find 0-flwe have 

d (<j> - Y') cos i - 1 



dt sin i 

and by proceeding exactly as above we shall find 



((o, sin + 1» 2 cos 0) 



«, cos + u>j sin = \ (A x + B z ) sin 2» - w 1 + J .4, - -Z? s sin w 1 



+ %A 2 -B l co&2n-n l +^A t + B, cosn 1 

sin (2» - 2» T ) 

also if « = i, + 2 where 2 denotes the sum of the quantities which occur 
in it, we shall have 

cos j - 1 cost- sin ».2- 1 cos t,- 1 -sin ^2 



sin <, + cos < 2 . / cos «, _, 

sin*, 1+- — - 2 
\ sin i, 



cost, -I / cos i,- 1 cosiA 

XJil_+ — sln^r") 



neglecting higher powers of 2 than the first, since these only are 
wanted. Multiplying together these factors, we shall have 



d (0 - ty) cos t, - 1 / —7— — =■ . r 

- K J± Y J = — r^ — iA l + B 2 sm2n-n i + t &c. 
dt suu, \ 

- (ATB2 A s - B t - A^By A^Bi) (l -, , s 

\ sin"*! 



cos*. - 1 COS I. 

+ ■ 



cos <! - 1 cos i,\ 2» 



sin *», ) 2n - n l 
Integrating and changing the signs on both sides 
cos i - 1 1 /A, + B 2 



sin 2m - 2re' 



<j)-e=—. — - ■' ; C os2w-» i +.. 

sin t, 2 \ 2» - m 1 



-\(A l + B,A 2 + B l -A 1 -B 2 A 2 -B l ) cos(2»-2n 1 ) 

» - n 1 2n - n 1 » 

We have now to substitute these values for t„ &c, in the function 

- - -. (1 + cos if sin 20- 26 
which occurs in the development of N. 
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In the first place, 



, _ 1 cos* (Aj+B 2 A t + B, - A,-B 2 A, -B. 

cos <= cos K + 2 + - - — — jt—7 

' 8 sin / \ ( 2n - n ) n 



sin In - 6m 1 



= cos(, {1 -42 s -sin«, 2-icos* — — 

1 T ' (2n-n 1 )n 1 



sin2w- 6» 1, 
Now, 

S* = - — — — - sin 2» - n 1 + - : — cos b, t 1 

2 2n-' 2 «' 



1 A, + B % A^B.-A,- B 2 A 2 - B, 
4 (2n-n 1 )n 1 



sin 2» - » l + t, 



P 



.: cos i = cos «, - sin i, 2 - x cos <, — tt— - sin 2» - 2m 1 

* (2n~n 1 )n 1 



where P = A, + U 2 , &c. 



1 + cos « = (1 + cos »,)* - 2 1 + cos <! sin <, 2 - 1 1 + cos <, cos i, 

P 

— — — . sin 2« - 2a 1 + sin < 2 ! 

(2n - n 1 ) n 

i> 



= 1 + cos«,*-21 + cosi sku,2 + (jsin s »-£ 1+cosi, cosi,) 
also 



(2»-n')n' 
sin ( 2n - n 1 ) 



sin20-20 = sin2»-2»' Li 008 /' l ( 2 2,v] 

( * 81111, ) 

+ 2 C ° 3 . t ' - cos 2n - 2w', ^i-i-P ■ cos' (2b- 2n') 

sin«, " 2n-n l n-n'n 

where 2, is the sum of the terms found in the first approximation. On 
expanding 22* , this becomes 



1 2» - 2» l |l— i^i — sm2n-2n\ +, &c. 

( sin« 2n-n l n 1 ) 

2cOSi-l -r— r „ /lcosi-1 1 1 



■ sin2»-2» 1 + — : cos 2»-2» 1 .2 1 - - — r 



sin t ' \4 sin *« (2» - n l ) n x 



1 n 

+ -. 



8 2» - n n - n l n l 
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If now we put for 2 . and 2! their values, and multiply this hy the 
value found above for 1 + cos . i J , we shall have, observing that the 

constant terms in 22, cos In - 'In 1 destroy each other, and retaining 

the constant term, 

,. . „j if 1 0O8«i- 1* COS li + l 1 

(1 +cos«)'sin2<£-20= ! r-r— ! +- 

(\ 4 8in s « 8 

. , 1= \ 1 1; , n \ 

sin \ - - 1 + COS »i cos «, 1 + COS H -r- — -- — rr-, P 

4 /2»-» 1 n 1 8 (2»-«)(«-»')» 1 ) 

the former part may be put into the form 

£ - ^ sin 'j - g 1 + cos « 2 cos ij 
or 

- | sin *< - | 4 cos t (1 + cos «, - 1 - cos »,) or - \ 1 + cos » 

therefore the whole expression becomes 

1 
(2» - n 1 ) (» - »') n 1 'In- «' »' 



- \ 1 + cos., 2 ( ta _ —z _n_. + sr-^r-il-P 



or 

- £ 1 + cos «,* - -— ; P. 

8 (» - »') n l 

Therefore, finally, the function which we have been examining contains 
the constant term 

frr) ? l\ TT ^ 12 oT^l? ^^ *** " ^ ir5 > 

or, as it may be more shortly written 

B-A\Z » l + cos, J 1 /. r ,^ B n 



C j 8 r" 4 (w - w l ) » : 

On putting for A x , &c, their values, and multiplying out, the latter 
factor becomes 



N(M y M\ . -, N 



N_ ( M,A-d , + M I U-£ 1 \ Ki 
A*B* 



Nl, r (C-B A-C\ „(A-G C-B 



+ ^{ M \ZkW + -WT)- M i-CB^-A^-)j nn - ni 

This may be much simplified for bodies nearly spherical, and might be 
put into the form 



NM y + M z ( , C-A C-B 



Ifi AB V ' AB 
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if we neglect small quantities which will be multiplied by B - A, &c), 
or observing that the latter factor is the same as D, it becomes simply 

N M X +M V 

lb AB ' 
and the term becomes 



B-A 3 1+cos*,'/. 1 N( w , r \ f\ 

which is true for every value of <. Suppose, for simplicity, that cos * 
does not much differ from unity ; then on this supposition we should 
have 

N ° "J ^X x ^ z - M > = T*\~8 W* 1 + *V^~ 2*^1 &C ' 
.-. M s +M y - £CjZ&i + H&x ~ 3 nV) = ~ (3 (*?«, + tfij - s, s ) 

upon substituting these values for 2V, and M x + J^, the expression a 
becomes 

2?- 4 135 / p Y 1 1 (3(^, + yA)-^') — f ., 

This term, it will be observed, is multiplied by the cube of ^, 
whereas the terms which have been previously examined were only 

multiplied by the square. If now we resume the expression in the de- 
velopment of N which is multiplied by 

r ( In - n 1 * v 2» - »') 

we shall see that, on continuing the approximation, {A x + jB 2 ) and 
(A % - Bi) will receive an increment such that when their new value is 
substituted for them in the above expression, it will be multiplied by 

it 
the cube of — . The term which so arises will not destroy that pre- 
viously found, but will be of the same order, and will modify it ; 
therefore it must be sought out. To find it, 
Let us return to the original equations for w> n> 2 , namely, 

du 1 3 p. C- B . (l+cos<- 

dt + ' = 2~r 3 ~~A~ Sm ' ( — 2 cosy-20-cos<cos0 

1 - cos ( ) 

r COS + 20/ 
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tf»2 


3 a A - C . ( 1 + COS » 

= 2r> B S1U ' 2 



sm 



<j>-ie-. ..} 



If, now, we take the terms already found in the first approximation 
for i, 0, vi}r, viz., 

1 A l -B 3 . At + B t 

« = 'i + 4 ; — sin» l + * : — oos n\ &c. 

and substitute them for i and 0, &c, in the above equations, we shall 
have terms whose argument is n - »'. And it will suffice for the pre- 
sent if we substitute the values in the first term which occurs on the 
right-hand side of each of these equations ; viz. that depending upon 
- 20 ; for, although that depending upon 0, that is the second, 
would produce terms of the same form, yet the increment so formed is 
best considered iu conjunction with the term in N which corresponds 
to it, namely, that depending upon Bin 20. Making, therefore, the sub- 
stitutions mentioned, we shall have to find the increment of u lt 

dw x 3 a G-B ,, , . , , lA t - B t . ; 

-57 + . . = - — — -. — (1 + oos«-4sin'») { ; sinm-w 1 

dt 8r, 4 a ( «' 

A, + B, 



+ ; — cosn-n' 



This might be integrated in the same way that the equations for 
«i, &c, were treated above ; but as it will suffice for the present pur- 
pose to reject quantities multiplied by products of 

G-B jC-A 
— and—, 

we may integrate it without introducing the second differential coeffi- 
cients. It will give 

3pG-B., . . ,_( A.-B, , 

"l" o ^ A~ (l + CO8l-j8Hl*0{- 7 -— COS»-»' 

A s + B, 



+ -. rr — ; sm n - w 

{n - n l ) »' 

which, if we make the same supposition as to » which was made before, 
becomes 



10, = — 



3 fL G-B I A t +B a A ^ £j 



cos n-n l + — — sin n - n* 



4r> A { n-n l n l {n-n')n l 

for the increment of u.\ above spoken of. 
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Similarly, for the increment of u>, -we have 

3 ft C-A ( A 3 + B, . A l -B 2 . , 

<o t = - — p {; - eos»-»' -— r— BUIB-H 1 

4 r* B ((» - n l ) n 1 (« - a 1 ) n l 

that is, A, is increased by 

3 /i C-B ( A x -B 2 



4 r 3 A \ (w - »') m 1 
and B % by 

3 fi C-A(A l -B 3 



4 r 5 B ^(rt-w 1 )^, 

whence .4, + -B s is increased by 

3 p. (C-A C-B\ A t -B, 



4 r 3 V S -4 7 (»-»') m 1 

where .4 i and -Bu on the right-hand side stand for the first approximate 
values of -4, and B 2 . 

Similarly A, - B z is increased by 

3 t [ °~ A - C ~^ \ A * + B * 
The increment, therefore, of N introduced by these, will be 



I ?15 3 (C-A C-B\{— 



A,-B 2 



Cr*~9 4\ B A j] iyil (2»?-« 1 ) (»-»»)n' 

1 — — A + 5, 



2 V l ;,1/ (2» -»,)(» -»,)»' 



which, if we put for A u &c, their first approximate values, and reject 
quantities depending upon higher powers of B - A becomes 



B-A 135 /> V 1 1 (3(xi\ +yi i z l )-z l 3 



-] *#!«!• (-B) 



J£C 32 \r»/ (2w-» 1 )ra 1 i^ 

It now remains to take account of the variations of <f> produced by 
those of «», 

# a. V 

(since ^ = * 3 , &c.) 

But with regard to these, it is easily seen that however they modify 
the preceding results, they cannot destroy them ; and for this simple 
reason, that every term in a', which will give rise to such termB will be 
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some quantity divided by C; and this, again, when substituted in the 
differential equations for to z will produce a quantity divided by C 1 , 
which therefore cannot destroy the terms above found. To proceed, 
then, we have in the first place, 

dw, B - A 3 a — s . . 

-± + &c.= c~8^ ' Sm ^~ )' 

which will give 

B-A S ft 2 . 

<0 S = » + — 7T- tj c . 1+cobi sin 2» - 2»' 

and hence 



<h=nt + — 7T- ^-1 + cos« sin 2»-2n 1 = n£ + P sin 2» - 2» l 
C7 r 3 



suppose, 



.-. cos = cosn + J P cos 3» - 2n' - ^ P cos » - 2»' 

multiply this by 

ib x = At cos « - n 1 , &e., 
and we have 



o, cos<ji= . . . + ^A X P cos2» -»' - £.4,Pcos 3»-2»' + &c. 



with terms of the same kind for « 2 Bin <f>, arising from B t sin n - n' 
.-. i»i cos <t> - o> sin = J ^4, - P, P (cos 2» - w 1 - cos 2» - Sre 1 ) 



. . _ „ / sin 2» - »' sm 3n - 2w' 

:l 1 + iJ 1 -P J P 



41 J V 2»-w J 2n-3» 1 



(the terms of the form -4, + B cos 2» - »', and may be rejected ; be- 
cause A t + P, contains P - 4 for a factor ; and as they are multiplied 
by P, whic h cont ains the same factor, they would produce terms mul- 
tiplied by B - A. Similarly, 



*-zr.±r A *-M-?zr-zr-* 



cos In - »' cos 2» - 3n l 



sin i * ' \ In - n l 2m - 3m 1 

with corresponding quantities for ijr. 

There is another way, also, in which terms wil l arise w hich must be 
taken into account, as follows : — The term P sin 2n - n 1 when intro- 
duced into the value of in the equations for «0[ and iu 2 will augment 
the value of A x so as to make it become A x (1 + ^P). Similarly, P 2 
will become P t (1 - 1 P), &c, and hence the term in 

A x + B. 



i i -zr r s iu 2» - »' 
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is augmented by 

A l -B l . 



4 tt r si 11 2» - n 1 ; 

and the terms in and ty are augmented in the same manner. Adding, 
then, these additional increments to those just found, we shall have for 
the whole increments depending upon P, 



, = , + A ( A-^ P cos (2» - ») + ( 4 i -^ ) -Pcos 2¥^3» l 
2 2W-W 1 v '* 2»-»' 

cos i , A, - B, „ . . cos < , A. — B. „ . , 

= »i = ^ l-J r 2 Psin2»-n' + -r— i-'— ^- Psin2»-3» I ) 

sin* a 2»-» 1 Bin « * 2n-3n' ' 



+ P sin 2» - 2k 1 

sm « * 2» - »' sin * * 2» - 3»' 

These terms must be added to those given in i, &c, by the first ap- 
proximation, and then substitute for *', &c, in N. Also it must be re- 
membered that corresponding to the term sin i sin 20 - 0, which arose 
from the multiplication together of sin i sin sin 20 - 20, there will 
be another term, sin < sin (20 - 80). Making the substitutions, there- 
fore, we shall have 

sin t (sin 20-5 - sin 20 - 30) = ( - + - 



\2 2« - »' 4 2n - 3»> 
1 



which will therefore give the term in N 

15/1 1_ 1 1 l\ 

2 \2 2n - n 1 + 4 2» - 3»> " n 1 ) l ~ Bl 

the last part of this, however, namely, that multiplied by will 

vanish ; for the term sin » sin sin (20- 20 l ), will produce in term in 
of the form QA 1 -B 1 cos 2» - n 1 , which, when introduced into the 
first part of JV, viz. 

B-A 



sin 20 - 20< 
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it is easily seen, will destroy the term spoken of. If, then, we cancel 
this term, and substitute for P, A, - B 2 , their values, we shall have, 

A 
observing that such a quantity as -^ may be replaced very nearly by 



unity, and rejecting any terms multiplied by products of B - A and 
C - A, &c, for this part of N 

n\n5S3JB-A 



(*) 



r 1 / 8 8 2 ABC \,2» - »' 2» - 3m 1 



h (»-«•) ) x ^i> — c 

The terms found are the only ones depending upon the same com- 
bination of a;/«„ &c, and having the same divisor D„ and not multi- 
plied by sin V„ or higher powers of B - A. 

There is one other term, however, i t will be advisable to take into 
account. It arises thus : the term P sin In - 2m 1 , when introduced into 
cos 8 — 6 l , &c, in the equations for *, and a> 2 , produces not only terms 
having n - n\ but also terms having 3» - 3»' for their arguments ; so 
that »„ &c, will contain terms of the form H cos Zn - 8» J , &c. Now, 
these terms will obviously contain the same combination of a?,**,, &c, 
in their coefficients, that A Jt &c, do; and for this reason they had 
better be retained. Their divisors, however, will be deficient ; but if 
we neglect quantities depending upon products of 

C-B C-A 
A B 

it is plain that the divisor of the latter may be pat equal to \ of the 
divisor of the former ; in fact, S will become A- A,P. 

In like manner if » 2 contains the term K sin Zn - 3n l , we shall have 
K= £ B 2 P. The effect of these terms will be to introduce into i the term 

4 J~ \ P sin 2i^Zn~\ 
12 2» - 3n l ' 

with similar terms in <j> and yfr ; and these, when introduced into the 
term - sin i sin 20 - 30 will give 

15 1 1 



PAt-B* 



2 12 2»-3n 1 

that is, they will reduce the term multiplied by 

1 
2n - 3»' 
in equation C to -Jrds the value it has there. 
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Suppose, now, we make this reduction, and collect the results, as 
given by equations (A t ), (i?,), and (C), and we have, finally, 

dw 3 _ 135 f /i\» B - A 1 1 / 1 1 \1 

dt = 32 [r> ) ABC D n - n l [y + 6.2«_3»V r*" 

It contains, therefore, a constant term ; and hence *> 3 will contain anew 
periodic term multiplied by the time. 

Corresponding also to equation (4), there will be a term arising 
from the second term in the first part of iV, viz., sin't sin 2p, and in 
(.B) a term arising from cos <p in the equation for «„ &c, as has been 
before mentioned : but these destroy each other. 

There will also be constant terms depending upon the arguments 

6\ 0-20', and 20 -C 1 , 

which arise in the same way as those given above ; that is, in every 
case the term arising from «,«, on the left side is identically de- 
stroyed by a term on the other side, while on the other hand the 
first term of the function iVwill produce, just as above, a term which 
is not destroyed, and another arising from the substitution of the 
next approximate values of*!, &c, in the second part ofiV; also in 
all cases the constant term produced by the function sin H sin 2<p in the 
first part of N is identically destroyed. The constant terms, however, 
in all these latter case3, are multiplied by sin f «, and depend upon quite 
different constant qualities. The one which has been examined ap- 
pears to be the most important, and is the only one which would exist 
independently of sia«; that is, when the plane of the equator does not 
differ sensibly from that of the orbit. The rest it is needless to say 
more about at present; and it only remains to discuss that already 
found. 

Discussion of the foregoing Results. 

(1.) Application to the Earth. 

The principal object there is to ascertain whether sufficient change 
has been produced in the motion of the earth to be perceptible by the 

observations of the last 2000 years. In the term -£■ which occurs, put 

r = a, and multiply numerator and denominator by /t + ^, where ^ is 
the mass of the earth, it becomes 

-j- ' or > P uttm s —z- = n • 



/* + /«, a- 
we have 



r 5 / XP + PiJ 



And for the moon, disturbed by the earth, we should have for the cop- 
responding quantity 
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n'° 



Also D " (n - n')\ nearly. 

(IbK 

The expression for —jf, may therefore be written 
dt 

B-A 135 / /t V n' 6 /l _J__W 

C 32 \/i + /. J (» - n l f \n l + 6 (2n - 3n') j a 8 

Where a, is the radius of the earth, a the radius of the orbit, to estimate 
the numerical value of this, we have 



U+*j U>) 



also — = —t ; as to the quantity — — 

a? 60 2 ^ * .4<fi 

it appears to be identical with one which occurs in the " Lunar Theory," 
and which has been estimated by M. Bessel from a great number of 
observations on the pendulum at about s^so ( see ^- Pontecoulant, 
" Systeme du Monde," vol. iv., p. 497). And we cannot, without mak- 
ing unwarrantable suppositions, attribute a greater value than this 

B-A 

either to — -~ — > or to the remaining quantity. Suppose, to fix the 
C 

ideas, that we were to put each of them at about ^j-jy. Also, the mul- 
tiplier 

n' 2 , »' n l n ri'n 

=- may be written — =-, or =- ; 

30 3 n 303 30 4 ' 

and to find the annual variations we must take n 1 to represent the mean 
motion in a year, i. e. about 17.000000", which, expressed in linear 
measure, may be put at about 85, so that we have for the annual 
variation of a 

135 1 85 _1_ 1 

"32* 80 1 30* (SO 1 3000 3 ' 

and the change in 2500 years will be thus multiplied by 2500, i e. about 
J_ 1 

70 (1000000) 3W * 

How, to produce an effiect such as is required to account for the accele- 
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ration of the moon, i. e. about \, it has been shown in the " Connais- 
sance du Temps," for 1800, that a variation in the length of the day of 
about a ten-millionth of a second, or, roughly, of about a million mil- 
lionth part of the entire length would suffice ; whereas, the value above 
given only amounts to about a seventy millionth part of this. And in 
like manner we shall find that the other terms mentioned, viz. those 
depending upon the argument l , will only produce quantities which 
may be put down as quantities of the same order, and apparently not 
more important than the above. The above effect is wholly due to the 
moon; that owing to the sun will be still smaller. Also, the earth 
has been treated as a solid body, no account being taken of its being 
partially covered with fluid. But Mr. Airy, I believe, has shown that 
the effect of the friction of the tidal wave in altering the rotation is 
quite insensible ; so that we may conclude that no effect sufficient to be 
perceptible in 2000 years has been produced by the action of gravity ; 
and hence direct perturbation of the rotatory motion is wholly inade- 
quate to explain the phenomena mentioned. 

A possible explanation, however, may be given by supposing that 
the slight alterations in figure which have been going on from geologi- 
cal or other causes may have been of a kind that shall suffice to cause 
such an alteration as would be required. If w is the angular velocity, 
it is well known that 



2mr 2 



e being a constant, m and r the mass and distance from the centre of any 
particle. Now, the mass remaining constant, it will be seen that a va- 
riation in length of the radius of gyration amounting to about a three 
thousandth part of an inch would be sufficient to produce the requisite 
effect. It is no violent supposition to suppose that such a change may 
have happened : perhaps the real difficulty would be, on the other 
hand, to suppose that a body with a radius of 4000 miles had remained 
so constant. But, on the other hand, the variation ought to be in in- 
crease, not a diminution in the radius ; and there is nothing to show 
that any changes that have taken place have tended in this direction 
rather than the other, further than this : — that any departure from the 
circular form in the sections parallel to the Equator tends to increase 
the radius of gyration ; so that if the general effect of geological changes 
has been to increase the inequalities of the surface, their effect would 
also have been slightly to lengthen the day, and hence, also, to cause 
an apparent acceleration of the moon. 

(2.) Application to the Moon. 

The chief question of interest which here occurs is : " what are the 
conditions of stable rotation?" or, in other words, what must be the re- 
lation between n and w 1 in order that there may be no permanent varia- 
tion in the rotation ? To answer this, it will be necessary to put down 

the general form which the equation for — will assume after the ap- 
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proximate values, including those arising from perturbation, have been 
substituted for <, <j>, and l in the general equation. "When this has been 
done, its general form will be 

da> 3 - 

-jt - N a + A x cos (» - n 1 1 + e) + B x sin (n - n l t + «) 



+ A 2 cos (2n - In t + 2e) + B 2 (sin In - 2n t + 2«) +, &c. 



+ P sin « M cos (p - q t + qe) + sin i t * * sin (p - q t - qe) +, &c. 

where N is the constant part, and where A„ &c, do not contain sin » 2 
as a multiplier, and where the last line represents the general form of 
those terms which are multiplied by sin i„ of which the terms in the 
expansion of N which have been used are an example, viz. 

/»' . 15 ._ , , . „ 

— sin » — - cos (2» - n l t + «), &c. 
f 3 2 

Let us omit for a moment the periodic terms, and consider the value 
of 2V , or rather that part of it found above, and which will be the only 
part when sin «, = 0. This term may be put under the form 

S 1 /l 1 \ 

+ 7T. 



Dn-n>\n x 6 2»-3»V 

Now this changes its sign, first, when the latter factor does, i. e. when 
12»- 7w, becomes ; that is, when n is something less than |n'. This, 
therefore, would, as far as thi3 term is concerned, be one condition 
under which the rotation would be stable ; but this would be no more 
than an approximate value, because it does not take account of terms 
multiplied by sin 2 «, which, for the moon, though small, is not 0. 
Again, it changes its sign when n - n 1 does; that is, supposing the co- 
efficient of to be positive, it would cause an acceleration or retar- 

» - n l 

dation according as n was greater or less than n 1 ; but in order that there 

might be stable rotation when n = »', it is evident that — 8 ought to 

change its sign by passing through zero, not by passing through infi- 
nity, as it appears at first sight to do when n - n 1 = 0. Let us examine 
what the true value is under such circumstances. It is quite evident 
that it must be either or infinite, since it is only by passing through 
one or other of these that it can change its sign. To see which it really 
is, it will be necessary to look back to the process by which the func- 
tion containing it was formed ; and if we do so, we shall see that the 
term is in reality only the first of a series, consisting of odd numbers of 
n - n\ having its signs alternately positive and negative ; so that 
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although, if we wer e to take only the first turn, it would certainly be 
infinite when n - n l = ; yet it does not follow that the entire series is 
so. On the other hand, it is quite evident that it is not ; for it arises 
from the multiplication together of such terms as sin t and sin n. Now, 
whatever be the value of « or x, sin t and x are never greater than unity, 
and hence their product must consist of terms whose sum cannot be 
greater than unity. In expanding sin x, &c, we only took the terms 
of lowest dimension that occurred, that is, we put 

» = •[ + o, or sin « = sin (i, + o) 

= sin « ]l - io* + r a* +j + cos » ( a - =- - &c. ) 

The higher powers of the quantities represented by a were rejected j if 
they had been retained, we should have had a series such as that men- 
tioned. If, then, the nonperiodic term cannot become infinite, and yet 
changes its sign when n becomes equal to «', it is plain that for such a 
value it must disappear, and that as far as this term only is concerned, 
n - n x = will be approximately a condition of stable rotation. If, 
however, we had taken into account terms depending upon the argu- 
ments 1 20 - 6', &c, it will be found that we should have introduced 
into the constant terms quantities multiplied by sin % and which do 
not change their sign when « = n 1 ; so that the entire value of N will 
be a quantity which does not change its sign, and cannot become when 
n=n x ; and the relation between n and n 1 thus obtained by equating 
iVJ, to will be the relation which ought to be used instead of n = »' ; 
but it would appear also that in addition to this relation, the conditions 
will also be very approximately satisfied by n = n\ To show this, let 
iV^.1 be the value which jV assumes when n = n 1 . 

do3% 
Then, the equation for — becomes 

dt 

-=~ = iV^,„i + A, cos e + 2?, sin e + A t cos 2e + B% sin 2e. 
A 2 may be rejected. Now, if we put this = it will give us a value by 

which € may be determined so as to satisfy the equation -^ = when 

at 
n = n' ; for, for no value of e can sin € and cos e be simultaneously equal 
to ; and since N^ 1 is very small compared with A t B u and espe- 
cially B 2 , it is evident that a possible value of e may be found to satisfy 
the equation 

Nton + A, cos e + &c. = 0. 

In other words, it must have a principal axis inclined at a particular 
angle to the radius vector of the disturbing body ; and this angle, 
though it appears to be small, cannot be 0. And these appear to be 
the only conditions when the mean value of » is nothing ; but in other 
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cases there will be other conditions, which will be seen thus : Let 

iV^) represent the value which iVi takes when the relation between 
n and n l is such that^m = qn 1 . Then, when this is the case, the equa- 
tion for a> 3 becomes 

fa 

"J! = -^(m) + -P sin i,""* sin (pe) + sin f cos (pe) + periodic terms. 

Now, the terms sin (pe) not containing the time, are constant quan- 

titles, and to find, when — iB 0, we must equate the three terms given 

above to zero, and this will give us an equation for determining what e 
must be in order that the body may revolve permanently with the relation 
pn = qn 1 . That such a condition may be possible, will, of course, imply 
that the coefficients, &c, in the above equation, have such a value as 
to give values of sin pe and cos pe not greater than unity. And it is 
evident, in order that such a requisition may be fulfilled, the coefficient 
of either of them cannot be small with respect to iV— . And this will 
show, that though there are several such relations, there can only be 
a limited number ; for, as the quantities p and q become large, so also 

does the power of — or — , which multiplies such terms, become 

large ; and hence for large values of p and q, the coefficients in ques- 
tion rapidly diminish, the more so when i is small, and there will not 
be many of them which are larger than Np^^, and consequently not 
many different conditions of stable rotation. Sow many there are it is, of 
course, impossible to say without more knowledge than we have, or ever 
can have, of the numerical values of the various quantities concerned. 

"We may conclude, then, certainly, that there will always be either 
an acceleration or retardation of a body revolving freely about a fixed 
point, and acted on by a/listurbing force moving round it, except when 
certain given relations exist between n and n 1 ; but which of the two it 
will be it is quite impossible to say, without knowing more about the 
form of the revolving body than we do of the moon ; for without such 
knowledge we cannot determine the algebraic signs of the various co- 
efficients. One thing, however, appears highly probable, and it is this : 
that if the conditions of equable rotation can be satisfied by values of 
» not very much greater than n 1 such as n = $ »', or, on the other hand, 
by values not very much less than n 1 , such as n = §»', if its rotation had 
ever been very much greater, or very much smaller than it is, it would 
seem that the change ought to have ceased when it came to a position 
of equable rotation, such as either of the former, without further 
diminishing or increasing till it became equal to w 1 ; and from hence 
it would appear that its rotation can never have been very different 
from what we actually observe it to be. 



